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ABSTRACT 

We present evolutionary calculations for very-low-mass stars and brown 
dwarfs based on synthetic spectra and non-grey atmosphere models which 
include dust formation and opacity, i.e. objects with T cff <2800 K. The 
interior of the most massive brown dwarfs is shown to develop a conductive 
core after ~ 2 Gyr which slows down their cooling. Comparison is made in 
optical and infrared color-magnitude diagrams with recent late-M and L-dwarf 
observations. The saturation in optical colors and the very red near-infrared 
colors of these objects are well explained by the onset of dust formation in the 
atmosphere. Comparison of the faintest presently observed L-dwarfs with these 
dusty evolutionary models suggests that dynamical processes such as turbulent 
diffusion and gravitational settling are taking place near the photosphere. As the 
effective temperature decreases below T e g m 1300 — 1400 K, the colors of these 
objects move to very blue near-infrared colors, a consequence of the ongoing 
methane absorption in the infrared. We suggest the possibility of a brown dwarf 
dearth in J, H, K color-magnitude diagrams around this temperature. 

Subject headings: stars: low-mass, brown dwarfs — stars: evolution 
— stars: colour - magnitude diagrams 
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1. Introduction 

Since the unambiguous identification of the first cool brown dwarf (BD) Gl 229B 
(Oppenheimer et al. 1995), the discovery of objects with temperatures cooler than the 
latest known M-dwarfs has been steadily increasing. These discoveries revealed two classes 
of objects requiring new spectral classifications. On one hand, a family of objects shows 
very red infrared (J-K, H-K) colors while the signature of metal oxides (TiO, VO), whose 
band strength index is used to classify M-dwarf spectral types, and hydrides (FeH, CaH 
bands) disappear gradually from the spectral distribution, as observed e.g. in GD 165B 
(Kirkpatrick et al. 1999a). Several of these so-called "L" dwarfs (Martin et al. 1999; 
Kirkpatrick et al. 1999b), have been discovered by the near-IR surveys DENIS (Delfosse et 
al. 1999) and 2MASS (Kirkpatrick et al. 1999b). On the other hand, a dozen of objects 
with properties similar to G1229B have been recently discovered with 2MASS (Burgasser et 
al. 1999), SDSS (Strauss et al. 1999) and the VLT (Cuby et al. 1999). The IR spectrum 
of these objects is characterized by the unambiguous signature of methane absorption 
in the H, K and L bands, the predicted dominant equilibrium form of carbon below a 
local temperature T « 1300 — 1500 K in the P « 5 — 10 bar range (Allard & Hauschildt 
1995; Tsuji et al. 1996; Fegley & Lodders 1996). This yields blue near-infrared colors, 
with J — K <0 but I — J>5 (see e.g. Allard 1999). These cool dwarfs are identified as 
"methane" dwarfs, sometimes called also "T-dwarfs" (Kirkpatrick et al. 1999b). 

The spectroscopic and photometric properties of L-dwarfs and, to less extent, of 
methane-dwarfs cannot be properly reproduced by dust-free atmosphere models, indicating 
that below T cS < 2800/T, theoretical models must include the formation and the opacity of 
dust grains (Tsuji et al. 1996, 1999; Ruiz et al. 1997; Jones and Tsuji 1997; Allard 1998; 
Leggett et al. 1998; Tinney et al. 1998; Basri et al. 1999; Kirkpatrick et al. 1999b). 

In previous papers, we have developed a consistent theory of the structure and the 
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evolution of low-mass, dense objects (Chabrier and Baraffe 1997, CB97) based on non-grey 
dust-free atmosphere models (Hauschildt et al. 1999, HAB99). Evolutionary models derived 
from this general theory describe successfully various observed properties of M-dwarfs 
down to the bottom of the main sequence: mass-magnitude relationships, color-magnitude 
diagrams, mass-spectral type relationships (Chabrier et al. 1996; Baraffe & Chabrier, 1996; 
Baraffe et al. 1997, 1998, BCAH98). The aim of the present paper is to extend these 
calculations into the brown dwarf domain by including the effect of grain formation both in 
the atmospheric equation of state (EOS) and in the opacity for objects with T eff < 2800K. 
The first comprehensive description of the effect of grain formation on the evolution of BDs 
was done by Lunine et al. (1989) and has been updated and incorporated in Burrows et al. 
(1993, 1997) calculations. However only a few grains were included in these calculations and 
the grain opacity is treated either as a frequency-independent Rosseland mean, implying 
grey atmosphere conditions, above 1300 K (see Figure 1 of Burrows et al. 1997), or is 
simply ignored (assuming grains settle immediately below the photosphere) below 1300 K, 
where non-grey atmosphere models are used. The present calculations extend significantly 
the number of grain species and include consistently their frequency-dependent opacity in 
the transfer equation, yielding consistent non-grey atmosphere structure, spectral colors 
and evolutionary calculations. The models are described in §2, whereas comparison between 
theory and observation in various color- magnitude diagrams (CMD) is presented in §3. The 
remaining uncertainties and shortcomings in the theory are discussed in the conclusion. 



2. Model description 



The main input physics involved in the present calculations has been presented in 
details in CB97 and BCAH98. The outer boundary conditions between the interior and the 
non-grey atmosphere profiles are described in CB97. Non-grey atmospheres are a necessary 
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condition when absorption coefficients strongly depend upon frequency, which is the case 
for all objects below T eff ~ 5000 K, the onset of molecular formation. The formation of 
grains in the present model atmospheres is described briefly in Leggett et al. (1998) and 
Allard (1999) and in details in Allard et al. (2000a). The equilibrium abundances of each 
grain species are determined from the Gibbs energies of formation, following the method 
developed by Grossman (1972). At each temperature and for each condensed phase under 
consideration, the variable K C (T) is calculated from the partial pressures of the species % 
forming the condensate (e.g. Al and O for corundum AI2O3), determined by the vapor 
phase equilibria, Pi = N{IZT, where iVj is the number of moles of species % and 1Z is the gas 
constant. This variable K c is compared to the equilibrium constant K eq , calculated from 
the Gibbs energy of formation of the condensate. The abundance of a condensed species 
is obtained by the condition that this species be in equilibrium with the surrounding gas 
phase, K c > K eq (Grossman, 1972). The opacities of the grains are calculated from Mie 
theory (see also Alexander & Ferguson, 1994). 

The formation of condensed species depletes the gas phase of a number of molecular 
species (e.g. VO, TiO, FeH, CaH, MgH) and of refractory elements such as Al, Ca, Ti, Fe, 
V (Fegley & Lodders, 1996; Allard 1998; Lodders 1999; Burrows & Sharp, 1999), modifying 
significantly the atmospheric structure and the emergent spectrum. For late M-dwarfs 
and for massive and/or young BDs, the main cloud formation is predicted to occur near 
the photosphere (see e.g. Allard 1998; Lodders 1999). The present calculations assume 
a grain-size distribution in the submicron range (Allard et al 2000a). This general grain 
treatment is consistent with the observations of the objects near the bottom of and below 
the MS, namely the DENIS and 2MASS objects (Delfosse et al. 1999; Kirkpatrick et al. 
1999b), GD 165B (Kirkpatrick et al. 1999a) and Kelu-1 (Ruiz et al. 1997), which all exhibit 
strong thermal heating and very red colors. Asides from this complex thermochemistry, 
grain formation is in fact a dynamical process and involves a balance between various 
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timescales, such as the condensation, evaporation, coagulation, coalescence and convection 
timescales (see e.g. Rossow 1978), not mentioning the fact that non-equilibrium species 
might be present, as for example in the atmosphere of Jupiter (Fegley & Lodders 1994). 
Since such dynamical processes are not incorporated in the theory at the present stage, we 
elected to conduct calculations under various extreme assumptions. 

The first set of models are based on the dust-free "NextGen" model atmospheres 
(Allard et al. 1996; HAB99) and are the same as in BCAH98. In the second set of models, 
all condensed equilibrium species are included both in the EOS and in the opacity, taking 
into account dust scattering and absorption in the radiative transfer equation. These 
models are referred to as "DUSTY" models. The third set of models includes the formation 
of grains in the atmosphere EOS, thus taking into account the photospheric depletion of 
dust-forming elements in the gas phase, but not in the transfer equation, therefore ignoring 
the opacity of these condensates. This case, referred to as "COND" models, mimics a rapid 
gravitational settling of all grains below the photosphere. This is similar to the Burrows et 
al. (1997) calculations below T e s < 1300 K (see above). 

As mentioned above, the DUSTY and COND models represent extreme situations 
which bracket the more likely intermediate case resulting from complex, and presently not 
understood, thermochemical and dynamical processes. 



2.1. Atmosphere profiles 

Figure [I] displays P-T atmosphere profiles for T e g = 1800 K, log g—5 and solar 
composition, under the afore- mentioned different assumptions for the dust treatment. The 
atmospheric heating due to the large grain opacity (the so-called greenhouse or backwarming 
effect) in the DUSTY model is clearly visible, and yields a significantly hotter structure. 
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The effect is dominant in the outer layers, where the spectrum forms, and decreases in 
the inner layers, where the connection with the interior profile occurs (see also Tsuji et al. 
1996). When dust opacity is not taken into account as in the COND models, the inner 
atmospheric structure is barely affected by the formation of dust and is similar to the 
dust-free case (NextGen) for this temperature. Only for significantly cooler T e g does dust 
formation start to affect substantially the inner structure compared to dust-free models. 

As indicated in Figure [l], convection occurs only in optically-thick regions (r > 1), 
in contrast to hotter objects (see BCAH98). However, this optical depth corresponds to 
only about one pressure scale height Hp, or even less, from the photosphere (r ~ 1) near 
which occurs the condensation temperature of most grains for objects near the bottom of 
the MS. Even for the cooler models (T c g ~ 1000 K) where convection retreats to deeper 
layers, the top of the convective zone lies only at a few Hp from the photosphere. This can 
have important consequences on the formation and settling of atmospheric grains. Indeed, 
turbulent diffusion, as produced for example by overshooting from the convection zone or 
by advection due to rotation, is a much more efficient mechanism than microscopic diffusion 
and sedimentation, at least for submicron-size particles. Although the temperature at the 
top of the convective zone is found to be generally above the condensation temperature 
of all grains, this turbulent diffusion will efficiently bring material upward to the region 
of condensation and maintain small-grain layers, which otherwise would have settled 
gravitationally, in this region. Note that for submicron-size species, the condensation time 
Tcond is much smaller than other dynamical timescales (see e.g. Rossow, 1978; Lunine et al. 
1989). As temperature decreases, the photospheric density increases, more species condense 
so that the particle density p p increases. Coalescence will become effective in growing large 
particles (r coa i/r con d oc T j p 2 p , Rossow 1978) and the afore-mentioned turbulent mixing could 
no longer prevent these large particles to fall out rapidly. Moreover the condensation line 
of some species (e.g. silicates) will drop below the photosphere (see e.g. Lodders 1999). 
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Within this general scheme, the DUSTY models correspond to a situation with very efficient 
turbulent mixing, whereas the COND models correspond to a situation where mixing is 
inefficient compared with sedimentation. Whether this general mechanism, and the related 
decrease of opacity, can explain the color saturation - and thus smaller backwarming - of 
the coolest L-dwarfs remains to be quantified correctly. Work in this direction is under 
progress. 

As anticipated from Figure [l], spectral and photometric properties are generally more 
affected by dust formation than the inner atmospheric structure and thus the evolution, as 
shown in the next section. This reflects the weak dependence of L and T c g upon opacity for 
BDs, with L oc ft ' 3 and T e g oc~ k oa (Burrows & Liebert, 1993). 

*** FIGURE 1 *** 

2.2. Evolutionary calculations 

Evolutionary calculations have been performed for 900 < T c g < 2800 K, which 
corresponds to masses between 0.1 M & and 0.01 M for ages 10 8 <t < 10 10 yr, following 
the method described in CB97. Figure 2 shows the evolution of T e g as a function of time 
for various masses and dust treatments. This illustrates the uncertainty expected from dust 
treatment on the cooling of BDs. For a given age, the differences between the two extreme 
DUSTY and COND models reach up to ~10% in T e g and ~25% in luminosity L. The 
hydrogen-burning minimum mass (HBMM) is only moderately affected by dust formation. 
The opacity due to the formation of different grains (e.g. silicates, Al 2 03, CaTiOs) produces 
a blanketing effect which lowers the effective temperature and luminosity for a given mass 
at the bottom of the main sequence, as noted initially by Lunine et al (1989). Models 
omitting grain opacity (NextGen and COND) yield m H BMM ~ 0.072 M Q , T e fj=1700K, 
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log L/Lq = 5 x 10~ 5 , whereas DUSTY models yield a slightly lower limit, ^h-hbmm ~ 0.07 
Mq, T e ff=1550 K, logL/L Q = 4 x 10 -5 , for solar composition. A lower luminosity requires 
a smaller H-burning nuclear energy in the core to reach thermal equilibrium, thus a lower 
central temperature and therefore a lower mass, which yields a lower HBMM. 

*** FIGURE 2 *** 

The present calculations incorporate new conductive opacities by Potekhin et al. 
(1999) in the interior. These opacities, initially developed for neutron star envelopes and 
white dwarf cores, improve previous calculations (Hubbard and Lampe, 1969; Itoh et 
al. 1983 and Mitake et al., 1984; Brassard and Fontaine, 1994) and cover the range of 
temperatures and densities characteristic of BD interiors. An interesting new property 
revealed by the present calculations is the growth of a conductive core for the most massive 
BDs. As these objects cool and contract, their interior becomes more and more degenerate 
(ip = kT/kTp 3.3 x 10~ 6 T (/i e /p) 2 / 3 , where Tp is the electron Fermi temperature and 
/i e is the electron mean molecular weight). Electron conductivity becomes more and more 
important and becomes eventually the main energy transport mechanism in the central 
regions, instead of convection. This occurs in the core of old (t > 2 Gyr) BDs in the mass 
range 0.03-0.07 M , with a maximum extension of the conductive core R con d/R-k ~ 0.6 — 0.7 
for the more massive ones, since the maximum central density occurs near the HBMM 
(Burrows et al. 1997; Chabrier and Baraffe 2000). In this mass range, the conductive core 
appears at t ~ 2-3 Gyr, for ijj <0.1, and increases rapidly, reaching up to 50%- 70% of the 
total mass at 10 Gyr for masses m « 0.05 — 0.06 M . Figure |3| shows the evolution of 
the conductive core Mamd/M* as a function of time for a 0.06 M & BD. We verified that 
these results still hold, although quantitatively different, when using previously published 
conductive opacities. These properties are essentially unaffected by the dust treatment in 
the atmosphere, since evolution depends weakly on it. 
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*** FIGURE 3 *** 

The onset and growth of a conductive core does affect the evolution. Once conduction 
sets in at the center, it increases dramatically the efficiency of energy transport compared 
to convection (v F » v conv , where v F = (2kT F /m e ) 1 ^ 2 is the electron Fermi velocity), 
decreasing the internal temperature gradient and yielding a cooler central temperature for a 
given energy flux T (VT 4 = Tf jR oc J-'.^cond, where n C ond is the conductive opacity). This 
increases electron degeneracy at the center (ip oc kT decreases). In first approximation, 
retaining only the zero-temperature electron contribution and the kinetic ionic contribution, 
the radius of a BD can be written R ~ R (l + tp) where R = 2.8 x 10 9 (-^)~ 1//3 /i~ 5/3 cm 
is the zero-temperature (fully degenerate, ip = 0) radius (Stevenson, 1991). Therefore the 
more degenerate conductive object contracts towards a slightly smaller radius, yielding a 
larger release of gravitational energy Sfl = J M P5(l/ p)dm within a Kelvin-Helmoltz time. 

This larger energy source increases the total binding energy of the star: 



where u is the specific internal energy, yielding a larger luminosity L for a given age, i.e. 
slowing down the cooling time r of the BD, since: 



where we have neglected the nuclear source of energy arising from some residual 
hydrogen-burning (see e.g. Figure 2 of Chabrier and Baraffe 2000). 

As an example, for a 10 Gyr old 0.06 M BD, T eff is 15% larger and L is 30% brighter 
when conduction is taken into account. An effect difficult to verify observationally, however, 




(1) 




(2) 
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given the faintness of such old BDs. Note that the onset of a conductive core occurs only for 
t >1 Gyr and is thus inconsequential for the application of the lithium-test, which concerns 
only younger objects (see e.g. CB97). 

3. Color - Magnitude diagrams 

In this section, we compare the present models with available observations of late-M, 
L- and methane-dwarfs in various CMDs. The comparison is limited to objects with known 
distances. Although trigonometric parallaxes have been determined for several 2MASS 
and DENIS L-dwarfs (Dahn et al. 2000), the only methane dwarf with known distance is 
G1229B. 

3.1. Optical CMD 

*** FIGURE 4 *** 

Figure |] compares models and various observations in a My-(V-I) CMD. Three 
L-dwarfs of the Dahn et al. (2000) sample are known to be very close binary systems 
(Martin, Brandner Sz Basri 1999; Koerner et al. 1999) with nearly equal mass components 
and their magnitude is corrected in the CMDs to take into account this property. Models 
with different dust treatments are displayed for ages of 0.1 Gyr (only for the DUSTY 
models; short-dash line) and 1 Gyr (solid line). For t > 1 Gyr, DUSTY isochrones are 
hardly distinguishable from one another in any of the following optical (VRI) and near-IR 
(JHK) CMDs, although the same masses correspond obviously to different magnitudes, 
and are not shown for sake of clarity. The properties of the DUSTY models for different ages 
are displayed in Tables 1-5. Note that we do not give the H-magnitude. Indeed, the present 
models overestimate the flux in this band and work remains to be done to understand 
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this shortcoming. We also give the lithium-abundance depletion for these models for 
applications of the lithium-test. The COND models displayed in Figure 4 (long-dash line) 
cover the mass range between 0.08 M Q (T eff = 2350K) and 0.03 M© (T eff = 1050^) at 1 
Gyr. 

For My brighter than 20 (T efr >2200 K), the NextGen models (dash-dot) are 
significantly bluer by 0.3-0.4 mag in (V — I) than the DUSTY models. This is due to 
the recent improved TiO line list (Schwenke, 1998) included in the DUSTY atmospheres 
(Allard et al. 2000a,b), while the NextGen models used the previous line list by Jorgensen 
(1994). As emphasized in BCAH98, the NextGen models predict too blue (up to ~ 0.5 
mag) optical colors compared to observations for 10 < My < 19 (see Figure 6 of BCAH98). 
As shown in Figure f|, the new TiO linelist improves significantly this situation (see also 
Allard et al. 2000b) and confirms the fact that the previous discrepancy stemmed primarily 
from a lacking source of opacity shortward of 1 /im (cf. §4.2 of BCAH98), and not from a 
major caveat in the treatment of convection. However, as shown in Figure |4] and Figure 
below, the situation is not completely satisfactory yet, with a remaining ~ 0.2-0.3 offset in 
colors at the bottom of the MS. 

*** FIGURE 5 *** 

The dust effect is clearly seen for My > 19, (V — I) >4.8, with a much more modest 
reddening as My increases (T e g decreases) for the DUSTY and COND models than for the 
dust-free models. This reflects metal and in particular TiO and VO depletion due to grain 
formation in the COND and DUSTY atmospheres. These species are two main sources 
of absorption in the optical in the temperature-range of interest. Note that, although the 
(V — I) colors are similar, the DUSTY models are significantly fainter than the COND 
ones for a given mass, revealing the grain opacity which is not included in COND models, 
and the depressed continuum in the optical due to the opacity-induced backwarming and 
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H 2 dissociation (Allard et al. 2000a). Interestingly enough, the same trend of change of 
slope in (V — I) colors at My ~ 19 appears in the available observations, but at a ~ 0.5 
mag bluer limit, i.e. (V — I) <4.5, than the theoretical predictions. This builds confidence 
in our basic treatment of grain formation but suggests that the gas depletion of the main 
absorbing species in the optical (TiO, VO, FeH) is underestimated in the present theory. 
Note in passing the limited interest of the V-band for BD search, with most massive BDs 
older than 1 Gyr being at My > 24. 

The saturation in optical colors is even more drastic in (R — I) (Figure |5|), where 
the DUSTY models yield almost constant (R - I) below T cff ~ 2300 K (Mj > 14) at 
(R — I) ~2.4-2.5. As seen in the figure, the distribution of various objects observed in 
these colors at the bottom of and below the main sequence, either in the field or in the 
Pleiades cluster, not only shows the same saturation effect as predicted by the models but 
is quantitatively in very good agreement with the predicted colors. As mentioned above, 
however, the ~ 0.2 — 0.3 mag offset above ~ 2300 K is still present. The COND models 
predict bluer colors as T e g decreases, and in general more flux in the V, R, I bands at a 
given T e g than the DUSTY models. This stems from (i) the absence of grain opacity in 
COND models and (ii) the backwarming effect in the DUSTY models (see Figure 1) which 
destroys IR-absorbing species such as e.g. H 2 , H 2 0, allowing more flux to escape at longer 
wavelengths. 



Although the L3.5 dwarf 2|mass-wj0035T591+182 (M v =21.15, M R =18.15, M I= 15.92 



i?-/=2.23)(Reid et al. 1999) seems to be lying on the COND track, it is probably misleading 
and reveals probable errors in the photometry. Indeed, the near-IR colors of this object are 
not reproduced by the COND models, but rather by the DUSTY ones (see Figure ^). 
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3.2. Infra-red CMD 

Figure |] displays the same comparisons in near-IR colors (J — K). As already 
mentioned by Leggett et al. (1998), the DUSTY models give the best fit to IR photometry 
of L- and very late M-dwarfs. They explain the very red colors characteristic of these 
objects in terms of the backwarming of the atmosphere due to grain absorption, yielding 
a reduction of the band strengths of the main IR absorbers, like e.g. H 2 0, as proposed 
initially by Tsuji et al. (1996). Such red colors are impossible to reproduce with dust-free 
opacity models (NextGen and COND). Note that for Mk brighter than ~ 11, the DUSTY 
models are slightly bluer in (J — K), by ~0. 1-0.2 mag, than the NextGen models. This 
arises from the use of different H2O line lists: Miller et al. (1994) in NextGen and the recent 
AMES list (Partridge and Schwenke, 1997) in DUSTY models. Although less complete for 
the higher energy transitions than the AMES line list, the Miller et al. (1994) list yields 
better agreement between models and IR photometric observations for T e s > 2200-2300 K 
(where dust does not affect the colors), since it relies on better potential surfaces for high 
temperatures (see Allard et al. 2000b for details). Below T e g- ~ 2000 K, these transitions are 
no longer important and the more complete AMES water linelist is to be used. Below this 
temperature, however, dust opacity is the main ingredient which shapes the IR colors, and 
the choice of the water line list is less crucial for photometric analysis. This illustrates again 
the remaining shortcomings in the present theory arising from still partially inaccurate 
molecular opacities, and stresses the need for future improvements in these inputs. 

*** FIGURE 6 *** 

When dust opacity is neglected (COND and NextGen), near-IR colors become bluer 
due to strong H2, H2O, CO and eventually CH4 absorption in the H and i^-bands. This 
yields very blue near-IR colors, with J — K^J — H <0.5, as illustrated by the good 
agreement of G1229B IR colors with the COND track (see also Burrows et al. 1997). The 
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effect is more pronounced for the COND models because of the lack of TiO absorption, 
which extends into the J-band at these very low effective temperatures. However, the 
COND models overestimate the flux at and shortward of ~ 1 \xm (V, R, I bands) for G1229B 
(Golimowski et al. 1998). The spectrum of GL229B indeed does not show strong indication 
for dust in its IR spectrum from 1 to 5 fim (Allard et al. 1996; Marley et al. 1996; Tsuji et 
al. 1996; Oppenheimer et al. 1998), but reveals discrepancies in the optical with dust-free 
models, as mentioned above. Several solutions have been proposed for this missing source of 
opacity in the optical, such as a specific haze opacity due to the presence of the nearby star 
(Griffith et al. 1998) or warm dust layers (Tsuji et al. 1999). These suggestions, however, 
are based on ad-hoc hypothesis required to fit G1229B spectrum and lack robust physical 
grounds. The methane-dwarf SDSS 1624 (Strauss et al. 1999), for example, does not have 
a hot companion, yet shows no flux blueward of 0.8 /im, indicating a strong absorption 
in the optical, as does G1229B. This suggests that the source of absorption is a common 
atmospheric BD property, precluding the haze opacity scenario. Strong alkali resonance 
lines (K I, Na I) have been proposed recently as an alternative, or complementary solution 
(Tsuji et al. 1999; Burrows et al. 1999). Clearly this question is not settled yet and requires 
futher observational tests and consistent calculations. 

For (J — K) ~ 1.0-1.9, good agreement is found between the data and the DUSTY 
models (see Figure §), except for the LO-dwarf 2 |mass-wj0746425| +200032 (Reid et al. 1999) 



(Mk=9.70, J-K=l.25), which lies significantly above both the other data and the models, 
making it a probable unresolved binary, as suggested by these authors. 

An interesting feature of this diagram is a hint that the redder L-dwarfs of Dahn et 
al. (2000), and the slope of the observed L-dwarf sequence itself, seem to deviate from the 
DUSTY tracks for Mk > 11.5 and (J — K) > 1.5, the reddest observed objects having a 
color (J — K) ~ 1.9, i.e. T e g ~ 1800 K. If the parallax and the photometry of these objects 
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are confirmed, this may very well illustrate gravitational settling of some grains below the 
photosphere for these objects, as mentioned in §2.1. This will yield less backwarming and 
thus progressively cooler atmospheres (see Figure 1), increasing H2O, CO and H2 absorption 
in the IR. Once CH4 forms near the photosphere, the integrated flux over the K-band is 
strongly reduced, although the flux still peaks around 1 jum, leading to G1229B-like colors. 
Whether this transition occurs sharply or smoothly remains to be elucidated. Note that 
the CO-CH4 transition occurs gradually with some of the two elements being present in the 
stability field of the other so that cool stars do contain some limited amount of methane 
and CO has been detected in G1229B (Noll et al. 1997; Oppenheimer et al. 1998). The 
recently discovered methane-dwarfs (Burgasser et al. 1999), unfortunately, do not bring 
information about this transition because of their color selection criterium, J — K < 0.4. 
In case of an abrupt change in color between "L" and "methane" -dwarfs, there should 
be indeed a "no-BD land" in the region 0.5 < J — K < 1.0 with a sharp transition for 
T e ff ~ T gas = 1300 — 1400 K, the temperature of formation of methane in the P ~ 3 — 10 bar 
range. This corresponds to masses m « 0.015 — 0.02 M for t=0.1 Gyr, m m 0.045 — 0.05 M 
for t=l Gyr and m ~ 0.065 M for t=5 Gyr. In any events, this concerns objects cooler 
than T e g w 1700-1800 K, the temperature derived with the DUSTY models along the 1 Gyr 
isochrone for LHS102B (Goldman et al. 1999; see also Figure 6), the reddest L-dwarf with 
known parallax, and for 2MASS-W1632+19, the faintest L-dwarf in the Dahn et al. (2000) 
sample with Mk ~ 12.75, (J — K) ~ 1.9. These temperatures agree within 100 K with the 
ones derived spectroscopically by Basri et al. (1999) based on alkali resonance absorption 
lines in the optical spectrum. Objects cooler than 2MASS-W1632+19 thus correspond to 
BDs with masses m<0.03 M for ages t > 0.1 Gyr, m<0.06 M for t > 1 Gyr and to 
essentially all BDs for t > 2.5 Gyr (see Figure 0). 

The same general behavior as in (J — K) is found in (J — H), with the same color 
saturation observed for the reddest L dwarfs for Mj > 13 and (J — H) < 1.2. Although not 
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shown in the present paper, a COND isochrone at t ~(a few)xl0 8 yr fits well GL229B both 
in (J — if), in contrast to Burrows et al. (1997) models, and (J — K) colors. Their models 
could not reproduce the L'-band of GL229B. Since then, Leggett et al. (1999) revised the 
flux of GL229B, with a significant modification of the V band, which is now reproduced 
by the Burrows et al. (1997) models. The reason of their poor fit in the H-band, however, 
compared with the present models, remains unclear. Indeed the explanation proposed in 
Burrows et al. (1997) was based on previous, less accurate observations of G1229B, as 
mentioned above, and does not hold anymore. 

Finally, Figure 7 displays L' versus (K — U). The K band is affected by H 2 , H 2 
and CH 4 absorption, whereas the V is strongly affected by methane (Allard et al. 1996; 
Marley et al. 1996; Tsuji et al. 1996). The trends displayed in the previous IR diagrams 
are replicated, with a severe reddening due to grain backwarming effects in the DUSTY 
models, but also a substantial reddening for the COND and grainless model, due to the 
strong CH 4 absorption in the K-band. 

The magnitudes of brown dwarfs and very-low-mass stars with dusty atmospheres 
for several isochrones are displayed in Tables 1-5. We stress, however, that the flux for 
these objects peaks around 1 fim, so that J,Z and H remain the most favorable bands for 
detection. 

*** FIGURE 7 *** 

4. Discussion and conclusion 

The comparisons in various CMDs displayed in the previous section illustrate the good 
general agreement between the present theory of evolution of very-low-mass stars and brown 
dwarfs with dusty atmospheres and observations. The models which include grain formation 
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and opacity in the atmosphere explain successfully the IR colors of late M- and L-dwarfs, 
as well as their blue loop in the optical, although improvement is still needed for the latter 
case to reach accurate quantitative agreement with observations. As shown by Kirkpatrick 
et al. (1999a) for GD165B and Ruiz et al. (1997) for Kelu-1, the DUSTY models improve 
also significantly the agreement with observed IR spectra. The present models provide a 
consistent description of the physics going on in the interior and the atmosphere of objects 
at the bottom of and below the main sequence, from M-dwarfs to L-dwarfs. They can be 
used with reasonable confidence in the IR bands (JKL'M) to calibrate the main properties, 
mass, age, T e ff, L of identified late-M and L-dwarfs. Shortcomings still appear in the optical, 
whatever the treatment of dust, in spite of noticeable improvement (see Figure 4). Baraffe et 
al. (1998) already pointed out a possible missing source of opacity in the optical to explain 
the observed (V — /) colors of M-dwarfs between T e R ~ 3600 K and ~ 2300 K. Although 
the new TiO linelist from AMES (Schwenke, 1998) improves appreciably the situation (see 
Figure |4j), a ~ 0.2 — 0.3 mag discrepancy still remains with observations (see Figure ||]). 
Below ~ 2300 K, comparison of the DUSTY models with observations in (V — I) (see 
Figure f|) strongly suggests either an underestimate of optical opacities or an underestimate 
of the depletion of molecular absorbers during grain condensation. Problems appear as well 
in spectroscopic analysis of atomic line profiles where the DUSTY models cannot match the 
optical observations (Tinney et al. 1998; Basri et al. 1999). Another remaining caveat of 
the present theory arises from shortcomings in the absorption coefficients of water, yielding 
a ~ 0.2 mag offset in (J — K) for T e g > 2300 K when using the most recent H2O line list 
(see e.g. Figure 6). Finally, the optical spectrum of GL229B remains to be described 
accurately. Present COND synthetic spectra overestimate the optical flux, pointing out a 
missing source of opacity. Work in this direction is under progress. 

Although it is important to solve these problems and eliminate these shortcomings, 
it is obvious from the various CMDs that substellar objects must be searched in the 
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near-infrared. Brown dwarfs around 1500 K radiate more than 90% of their energy at 
wavelengths longward of 1 /im. The flux peaks at 1.1 and 1.3 jum, and J,Z,H are the 
favored broadband filters for detection. It is thus reinsuring, given the complexity of 
the underlying physics, that the present models of dusty substellar objects describe the 
observations in near-IR bands with reasonable success. 

If dust formation and absorption successfully explain the transition between M- and 
L-dwarfs, a theoretical and observational gap still exists between L- and methane-dwarfs. 
Models including dust opacity reach rapidly very red IR colors (J — K) > 2, whereas 
observations of the reddest or faintest L dwarfs show a saturation of their (J — K) and 
(J — H) colors at 1.9 and 1.2 respectively, which corresponds to T e g ~ 1700-1800 K (see 
Figure |6|). This very likely illustrates ongoing dynamical processes near the photosphere. 
As outlined in §2.1, the proximity of the top of the convection zone from the photosphere 
may induce turbulent mixing and prevent the grains to diffuse downward. As convection 
retreats to deeper layers for cooler objects, gravitational settling of condensed species may 
lead to cooler atmospheres, because of the decreasing opacity, with a transition to very blue 
near infrared colors once methane forms near the photosphere. These dynamical processes 
of grain formation and diffusion are not included in the theory yet and represent indeed the 
next major challenge. Theory definitely needs some guidance from observations through 
this puzzle, by getting spectra, parallaxes and photometry of the coolest L-dwarfs and by 
either filling or confirming the presently expected brown dwarf dearth between the coolest 
L-dwarfs and the hottest methane dwarfs around T e g « 1300 — 1400 K. 

Note: Various isochrones for the DUSTY models from 10 7 yr to 10 10 yr, covering 
a range of mass from 0.01 M & to 0.1 M Q and T e g from 900K to 3000K are available by 
anonymous ftp (same format as Table 1): 

ftp ftp.ens-lyon.fr 
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username: anonymous 

ftp > cd /pub/users/ CRAL/ibaraffe 

ftp > get DUSTY00_models 

ftp > quit 

Note that the DUSTY models are given in the Tables down to T c fj=900K, but this 
limit is clearly unrealistic for dusty atmospheres, as seen in the various figures of the paper. 

Special requests for models in any particular filters can be addressed to I. Baraffe. 
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Table 1: DUSTY model isochrone for 0.1 Gyr. T e ff is in K, the gravity g in cgs. The lithium 
abundance is normalized to the initial abundance, which is Xu = 10~ 9 by mass fraction in 
the present calculations. The VRI magnitudes are in the Johnson-Cousins system (Bessell 
1990), JHK in the CIT system (Leggett 1992), L' in the Johnson-Glass system and M in the 
Johnson system. 



m/M Q 


T ef.f 


log L/Lq 


log g 


Rj R Q 


Li/Lio 


My 


Mr 


Mi 


Mj 


M K 




M M 


0.010 


916. 


-4.98 


4.21 


0.130 


1.0000 


42.64 


35.06 


31.45 


22.06 


15.20 


11.71 


11.30 


0.012 


1333. 


-4.25 


4.22 


0.141 


1.0000 


29.33 


24.92 


22.13 


16.50 


12.14 


10.10 


10.12 


0.015 


1290. 


-4.37 


4.37 


0.132 


1.0000 


30.63 


25.96 


23.08 


17.14 


12.52 


10.32 


10.36 


0.020 


1441. 


-4.21 


4.53 


0.127 


1.0000 


27.30 


23.35 


20.73 


15.68 


11.86 


10.14 


10.20 


0.030 


1878. 


-3.75 


4.71 


0.127 


1.0000 


21.21 


18.37 


15.96 


12.22 


10.76 


9.74 


9.89 


0.040 


2234. 


-3.41 


4.80 


0.132 


1.0000 


18.72 


16.41 


13.95 


11.04 


10.16 


9.37 


9.62 


0.050 


2483. 


-3.18 


4.84 


0.140 


0.9990 


17.07 


15.14 


12.86 


10.47 


9.70 


9.04 


9.28 


0.055 


2576. 


-3.09 


4.86 


0.145 


0.9920 


16.44 


14.65 


12.48 


10.27 


9.51 


8.90 


9.12 


0.060 


2652. 


-3.01 


4.87 


0.149 


0.9630 


15.95 


14.26 


12.18 


10.11 


9.36 


8.77 


8.98 


0.070* 


2768. 


-2.88 


4.88 


0.159 


0.6680 


15.19 


13.64 


11.70 


9.83 


9.09 


8.55 


8.72 


0.072 


2787. 


-2.86 


4.88 


0.161 


0.5380 


15.07 


13.54 


11.63 


9.78 


9.04 


8.51 


8.68 


0.075 


2814. 


-2.83 


4.88 


0.164 


0.3290 


14.89 


13.40 


11.52 


9.71 


8.97 


8.45 


8.62 


0.080 


2854. 


-2.78 


4.89 


0.169 


0.0664 


14.65 


13.19 


11.36 


9.61 


8.87 


8.36 


8.52 


0.090 


2920. 


-2.69 


4.89 


0.178 


0.0000 


14.23 


12.84 


11.09 


9.41 


8.68 


8.20 


8.34 


0.100 


2971. 


-2.63 


4.90 


0.185 


0.0000 


13.92 


12.58 


10.89 


9.27 


8.54 


8.07 


8.21 



* HBMM 
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Table 2: 


Same 


as 1 


rable 


1 for 0.5 


Gyr 


















m/M Q 


T ef.f 


loj 


lL/L 


log 9 


R/R Q 


Li/Li 


My 


Mr 


Mi 


Mj 


M K 


M L , 


M M 


0.020 


901. 




-5.14 


4.65 


0.111 


1.0000 


43.80 


35.89 


32.23 


22.56 


15.59 


12.11 


11.67 


0.030 


1170. 




-4.73 


4.86 


0.106 


1.0000 


33.93 


28.57 


25.49 


18.78 


13.65 


11.09 


11.10 


0.040 


1461. 




-4.36 


5.01 


0.104 


1.0000 


27.41 


23.50 


20.90 


15.94 


12.20 


10.51 


10.62 


0.050 


1751. 




-4.06 


5.11 


0.103 


0.9540 


22.45 


19.47 


17.10 


13.18 


11.35 


10.23 


10.36 


0.055 


1902. 




-3.91 


5.15 


0.103 


0.5610 


21.42 


18.64 


16.25 


12.60 


11.20 


10.16 


10.32 


0.060 


2048. 




-3.78 


5.19 


0.104 


0.0069 


20.56 


17.91 


15.46 


12.12 


11.00 


10.04 


10.25 


0.070* 


2295. 




-3.56 


5.23 


0.106 


0.0000 


18.95 


16.62 


14.18 


11.42 


10.58 


9.78 


10.06 


0.072 


2332. 




-3.52 


5.24 


0.107 


0.0000 


18.71 


16.44 


14.01 


11.34 


10.51 


9.74 


10.02 


0.075 


2400. 




-3.46 


5.24 


0.108 


0.0000 


18.26 


16.09 


13.71 


11.18 


10.39 


9.65 


9.93 


0.080 


2505. 




-3.37 


5.25 


0.111 


0.0000 


17.56 


15.56 


13.27 


10.95 


10.18 


9.51 


9.79 


0.090 


2678. 




-3.19 


5.24 


0.118 


0.0000 


16.40 


14.67 


12.58 


10.57 


9.83 


9.24 


9.48 


0.100 


2809. 




-3.05 


5.23 


0.128 


0.0000 


15.52 


13.97 


12.06 


10.25 


9.52 


8.99 


9.19 



* HBMM 
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Table 3: Same as Table 1 for 1 Gyr 



ui/Mq 


T ef.f 


log L/Lq 


log g 


Rj Rq 


Li/Lio 


My 


Mr 


Mi 


Mj 


M K 




M M 


0.030 


958. 


-5.11 


4.90 


0.102 


1.0000 


40.76 


33.66 


30.17 


21.50 


15.22 


11.98 


11.68 


0.040 


1189. 


-4.77 


5.06 


0.098 


1.0000 


33.53 


28.29 


25.25 


18.69 


13.69 


11.20 


11.24 


0.050 


1424. 


-4.48 


5.18 


0.096 


0.9320 


28.26 


24.20 


21.54 


16.42 


12.53 


10.74 


10.86 


0.055 


1550. 


-4.34 


5.22 


0.095 


0.3790 


26.08 


22.48 


20.00 


15.39 


12.08 


10.59 


10.71 


0.060 


1675. 


-4.21 


5.27 


0.094 


0.0000 


23.95 


20.72 


18.33 


14.14 


11.71 


10.48 


10.59 


0.070* 


2012. 


-3.88 


5.32 


0.096 


0.0000 


21.10 


18.37 


15.95 


12.47 


11.24 


10.24 


10.43 


0.072 


2089. 


-3.80 


5.32 


0.097 


0.0000 


20.51 


17.87 


15.41 


12.17 


11.10 


10.15 


10.37 


0.075 


2185. 


-3.71 


5.33 


0.098 


0.0000 


19.82 


17.32 


14.85 


11.85 


10.92 


10.03 


10.29 


0.080 


2358. 


-3.54 


5.32 


0.103 


0.0000 


18.66 


16.41 


14.00 


11.38 


10.57 


9.80 


10.08 


0.090 


2625. 


-3.27 


5.28 


0.113 


0.0000 


16.81 


14.99 


12.84 


10.74 


9.99 


9.38 


9.64 


0.100 


2790. 


-3.08 


5.25 


0.125 


0.0000 


15.68 


14.10 


12.16 


10.33 


9.59 


9.05 


9.26 
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Table 4: Same as Table 1 for 5 Gyr 



itl/Mq 


T ef.f 


log L/Lq 


log g 


Rj Rq 


Li/Lio 


My 


Mr 


Mi 


Mj 


M K 


M v 


M M 


0.050 


907. 


-5.35 


5.26 


0.086 


0.9210 


43.03 


35.38 


31.77 


22.53 


15.95 


12.57 


12.11 


0.055 


990. 


-5.21 


5.32 


0.085 


0.2840 


39.43 


32.76 


29.36 


21.21 


15.26 


12.23 


11.95 


0.060 


1097. 


-5.04 


5.37 


0.084 


0.0000 


35.60 


29.92 


26.76 


19.74 


14.50 


11.81 


11.74 


0.070* 


1593. 


-4.36 


5.40 


0.087 


0.0000 


25.54 


22.08 


19.67 


15.23 


12.13 


10.73 


10.84 


0.072 


1754. 


-4.18 


5.40 


0.089 


0.0000 


22.96 


19.93 


17.58 


13.63 


11.68 


10.54 


10.66 


0.075 


1998. 


-3.92 


5.38 


0.092 


0.0000 


21.31 


18.55 


16.14 


12.62 


11.35 


10.32 


10.50 


0.080 


2295. 


-3.61 


5.34 


0.100 


0.0000 


19.12 


16.77 


14.32 


11.57 


10.73 


9.91 


10.20 


0.090 


2621. 


-3.28 


5.29 


0.113 


0.0000 


16.85 


15.02 


12.86 


10.75 


10.01 


9.39 


9.65 


0.100 


2791. 


-3.08 


5.25 


0.124 


0.0000 


15.68 


14.10 


12.16 


10.33 


9.59 


9.05 


9.26 



* HBMM 



Table 5: 


Same 


as Table 1 for 10 Gyr 


















m/M Q 


T ef.f 


log L/L© 


log g 


Rj Rq 


Li/Lio 


My 


Mr 


Mi 


Mj 


M K 


M L , 


M M 


0.060 


971. 


-5.29 


5.40 


0.081 


0.0000 


40.34 


33.44 


30.00 


21.61 


15.52 


12.43 


12.09 


0.070* 


1552. 
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FIGURE CAPTIONS 



Fig. 1. — Atmosphere profiles for T eff =1800 K and log g—5. Different treatments of dust 
are considered. Solid line (full circle): DUSTY; long-dashed line (full square): COND; dash- 
dotted line (open circle): NextGen (see text). The photosphere (r Rossc i and ~ 1) is located 
at T ~ 2100 K for all models, whereas the end of the curves correspond to r = 100, the 
boundary condition with the interior profile. The symbols on each curve indicate the onset 
of convection 

Fig. 2. — Evolution of T e ff as a function of time (in yr) for different masses and atmosphere 
models. Solid lines: DUSTY; long-dashed lines: COND; dash-dotted lines: NextGen. Masses 
are indicated in M Q 

Fig. 3. — Evolution of the conductive core M cond /M tot (shaded area) as a function of time 
for a 0.06 M brown dwarf. 

Fig. 4. — CMD for the optical colors (V — I) at the bottom of the Main Sequence and below. 
The data is from Monet et al. (1992) (filled circles), Dahn et al. (1995) (crosses), Dahn 
et al. (2000) (filled squares) and Reid et al. (1999) (filled triangles). Kelu-1 (Ruiz et al. 
1997) is indicated by a full diamond. Theoretical models correspond to various isochrones: 
short-dashed and solid lines: DUSTY models for 0.1 and 1 Gyr respectively; long-dashed 
line: COND models for 1 Gyr; dash-dotted line: NextGen models for 1 Gyr. Masses (in 
Mq) and T cff are indicated on the 1 Gyr DUSTY isochrone by open circles. The COND 
models go from 0.08 M (same as DUSTY) to 0.03 M . See Tables 1-5 for the location of 
the HBMM at different ages 
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Fig. 5. — CMD for the optical color (R — I). The filled circles are the Pleiades data from 
Bouvier et al. (1998). The other symbols and the models are the same as in Figure See 
Tables 1-5 for the location of the HBMM at different ages 



Fig. 6. — CMD for the near-IR color (J—K). The dots are the data from Leggett (1992). The 
other symbols are the same as in Figure [|. Some identified BDs are indicated by diamonds: 
Kelu-1 (Ruiz et al. 1997); GD165B (Kirkpatrick et al. 1999a); LHS102B (Goldman et al. 
1999). Theoretical models are the same as in Fig |4]. Crosses on the COND track indicate 
masses in M Q and T e g- in K. See Tables 1-5 for the location of the HBMM at different ages 



Fig. 7.— CMD for the IR color (K — L'). The dots are data from Leggett (1992). Theoretical 
models are the same as in Figure |6|, with the same masses indicated on the COND isochrone 
(long-dash line, x). The dotted curve corresponds to the DUSTY isochrone for 10 Gyr. 
Masses in M & are indicated on the DUSTY isochrone for 1 Gyr (solid line, O); T c g for 
m > O.O8M are the same as in Figure |j. See Tables 1-5 for the location of the HBMM at 
different ages 
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